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Azimuthal asymmetries in high-energy processes, most pronounced showing up in combi- 
nation with single or double (transverse) spin asymmetries, can be understood with the 
help of transverse momentum dependent (TMD) parton distribution and fragmentation 
functions. These appear in correlators containing expectation values of quark and gluon 
operators. TMDs allow access to new operators as compared to collinear (transverse mo- 
mentum integrated) correlators. These operators include nontrivial process dependent 
Wilson lines breaking universality for TMDs. Making an angular decomposition in the 
azimuthal angle, we define a set of universal TMDs of definite rank, which appear with 
process dependent gluonic pole factors in a way similar to the sign of T-odd parton dis- 
tribution functions in deep inelastic scattering or the Drell-Yan process. In particular, 
we show that for a spin 1/2 quark target there are three pretzelocity functions. 

Keywords: Parton distributions; Transverse Momentum Dependence; QCD. 

PACS numbers: 12.38.-t, 13.85.Ni, 13.85.Qk 



1. Introduction 

To study the connection between partons and hadrons in high energy processes 
through parton distribution functions (PDF) and parton fragmentation functions 
(PFF), the starting points are forward matrix elements of parton fields, such as the 
quark-quark correlator 

*«(p|p) = / (01 e ^ (p\^(o)MO\p), (i) 

where a summation over color indices is understood. This replaces the correlator 
$ oc (j> + m) for a single incoming fermion. In the case of hadrons one also needs 
quark-quark-gluon correlators such as 

^ Mj (p-pun\p) = J e i(p " PlH e*™ (P^A^MOIP)- (2) 



QCD Evolution Workshop 2012 Buffing-Mulders 



2 M. G. A. Buffing and P. J. Mulders 

The basic idea is to factorize these hadronic (soft) parts in a full diagrammatic 
approach and parametrize them in terms of PDFs. This requires high energies in 
which case the momenta of different hadrons obey P-P' oc Q 2 , where Q 2 is the hard 
scale in the process. In that case the hadronic momenta can be treated as light-like 
vectors P and the hard process brings in a conjugate light-like vector n such that 
P-n = 1, for instance n = P 1 /P-P 1 . One makes a Sudakov expansion of the parton 
momenta, 

p = xP + p T + (p-P - xM 2 )n, (3) 

with x = p + = p-n. In any contraction with vectors outside the correlator, the 
component xP contributes at order Q, the transverse component at order M ~ Q° 
and the remaining component contributes at order M 2 /Q. This allows consecutive 
integration of the components to obtain from the fully unintegrated correlator in 
Eq. ((TJ the TMD light-front (LF) correlator 



the collinear light-cone (LC) correlator 



£-ra=0 



(4) 



(5) 



£-n=£ r =0 Or £ 2 =0 



or the local matrix element 

4> 



y- TO J (O)^(0|i > >| e=o . (6) 



The importance of integrating at least the light-cone (minus) component p~ = p-P 
is that the expression is at equal time, i.e. time-ordering is not relevant anymore for 
TMD or collinear PDFs - . For local matrix elements one can calculate the anoma- 
lous dimensions, which show up as the Mellin moments of the splitting functions 
that govern the scaling behavior of the collinear correlator $(a;). We note that the 
collinear correlator is not simply an integrated TMD. The dependence on an up- 
per limit <fr(x;Q 2 ) = J d 2 p T <&(x,p T ) is found from the anomalous dimensions 
(splitting functions). One has an a s /p% behavior of TMDs that is calculable using 
collinear TMDs and which matches to the intrinsic non-perturbative p T -behavior^l 
We note that in operator product expansion language, the collinear correlators in- 
volve operators of definite twist, while TMD correlators involve operators of various 
twist. 



2. Color gauge invar iance 

In order to determine the importance of a particular correlator in a hard process, 
one can do a dimensional analysis to find out when they contribute in an expansion 
in the inverse hard scale. Dominant are the ones with lowest canonical dimension ob- 
tained by maximizing contractions with n, for instance for quark or gluon fields the 



QCD Evolution Workshop 2012 Buffing-Mulders 



Generalized Universality for TMD Distribution Functions 3 

minimal canonical dimensions dim[-0(O)^. ip(£)] = dim[F na (0) -F n/3 (£)] = 2, while an 
example for a multi-parton combination gives dim[?/>(0)?4 -A™ (77) "0(C)] = 3. Equiva- 
lently, one can maximize the number of P's in the parametrization of Of course 
one immediately sees that any number of collinear n-A(rj) — A n (rj) fields doesn't 
matter. Furthermore one must take care of color gauge invariance, for instance 
when dealing with the gluon fields and one must include derivatives in color gauge 
invariant combinations. With dimension zero there is iD n = id 71 + gA n and with 
dimension one there is iD" = id" + gA". The color gauge-invariant expressions for 
quark and gluon distribution functions actually include gauge link operators, 




U m =Vexp[-i I dC^(0], (7) 
connecting the nonlocal fields, 



x Tr (P,S\ F n »(0) U m F™(0 U{ m \P,S) 



(8) 



LF 



(9) 



LF 



For transverse separations, the gauge links involve paths running along the mi- 
nus direction to ±00 (dimensionally preferred), which are closed with one or 
more transverse pieces at light-cone infinity. The two simplest possibilities are 
t/M = E7|q ±00 , E7jQ r £ i ^[±oo i leading to gauge link dependent quark TMDs 

i'q 4 ' (x,p T ) I^El. For gluons, the correlator involves color gauge-invariant traces of 
field-operators F na , which are written in the color-triplet representation, requiring 
the inclusion of two gauge links £/[o,£] and UL j . Again the simplest possibilities are 
the past- and future-pointing gauge links giving even in the simplest case four 
gluon TMDs T^' ±] \x,p T ). 

Using the dimensional analysis to collect the leading contributions in an expan- 
sion in the inverse hard scale, one will need the above quark and gluon TMDs for 
the description of azimuthal dependence. Taking the Drell-Yan (DY) process as an 
example, one can look at the cross section depending on the (small!) transverse 
momentum q T of the produced lepton pair, 

a(xi,x 2 ,q T ) = J d 2 p lT d 2 p 2T S 2 (p lT + p 2T - q T ) 

X §Y ] {xi,Pi t )$ 2 ^\x2,P2t)v(x 1 ,X2,Q), (10) 



which involves a convolution of TMDs. What is more important, it is the color flow 
in the process, in this case neutralized in the initial state, that determines the path 
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in the gauge link in the TMDs, in this case past-pointing ones. In contrast in semi- 
inclusive deep inelastic scattering one finds that the relevant TMD is with a 
future-pointing gauge link. In a general process one can find more complex gauge 
links including besides Wilson line elements also Wilson loops. In particular when 
the transverse momentum of more than one hadron is involved, such as e.g. in the 
DY case above, it may be i mpos sible to have just a single TMD for a given hadron 
because color gets entangled 

T he co rrelators including a gauge link can be parametrized in terms of TMD 
PDF S depending on x and p T , 

{PtSt 
jf _ f^ U \x,p T ) fZ_ +gM( X ,p T ) l5 

+ h ^ {x y T ) l a T + hi} U \ X ,p T )^ + i h^ U \ X ,p T )^^, (11) 

with the spin vector parametrized as S 7 ' = S L P^ + St + M 2 S L n^ and shorthand 
notations for gf^ and hf} U ^ , 

= S L gM(x,pl) - E^gM^pl). (12) 

For quarks, these include not only the functions that survive upon p T -integration, 
fi(x) = q{x), g\{x) = Aq(x) and h\ (x) — 5q(x), which are the well-known collinear 
spin-spin densities (involving quark and nucleon spin) but also momentum-spin den- 
sities such as the Sivers function f 1 ^(x,p T ) (unpolarized quarks in a transversely 
polarized nucleon) and spin-spin- momentum densities such as g^j<{x,p T ) (longitu- 
dinally polarized quarks in a transversely polarized nucleon). 

The parametrization for gluons, following the naming convention in Ref. [9| is 
given by 

PtSt 

2xT^(x,p T ) = ff\x,pl) + «^^_ f^ 9[U] (x,pl) 

%^f- hff u \x,p T ) - " T bT ^ T Pt h°V\x,pl). ( 



3. Transverse moments 

In many cases, it is convenient to construct moments of TMDs in the same way as 
one considers moments of collinear functions. For §(x) in Eq. ((5| one constructs 
moments 



2n 



LC 

(14) 

LC 
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Integrating over x one finds the connection of the Mellin moments of PDFs with 
local matrix elements with specific anomalous dimensions, which via an inverse 
Mellin transform define the splitting functions. Similarly one can consider transverse 
moment weighting starting with the light- front TMD in Eq. (HJ), 

(2tt)* 



x (p\m t/[™ ±00] u[ 0tXt] iD«(±oc) u^mip) 



.(15) 

LF 

Integrating over p T gives the lowest transverse moment, which appears in the q T - 
weighted result of Eq. (|10[) . This moment involves twist-3 (or higher) collinear multi- 
parton correlators, in particular the quark-quark-gluon correlator 



{2-Kf 

x e^o (P\m U$ tV] F na (rj) U^ (] ^)\P) 



• (16) 

LC 



In terms of this correlator and the similarly defined correlator <& D ( x ~ x i ■> x i \ x ) one 
finds 



d 2 p T p«$M(x,p T ) = $g(x) +4 C/I 7r$g(x), (17) 



with 



= J dxi $n(x — xi,xt\x) — J dxiPV— §p a (x - xi,xi\x), 
= <f> F a (x,0\x). 

The latter is referred to as a gluonic pole or ETQS-matrix element They 
are multiplied with gluonic pole factors C^ 7 ' (e.g. Cq = ±1), that tell us that 
new functions are involved with characteristic process dependent behavior l l^ l l'H 
This behavior is for the single transverse moments also coupled to the behavior 
under time-reversal. While $g is T-even, $g is T-odd. Since time-reversal is a good 
symmetry of QCD, the appearance of T-even or T-odd functions in the parametriza- 
tion of the correlators is linked to specific observables with this same character. In 
particular single spin asymmetries are T-odd observables. 

The weighting with transverse momenta can also be analyzed by studying the 
parametrization in PDFs. For single p T -weighting, only PDFs with one prefactor of 
p T in the parametrization in Eq. (jllj) survive. The <&q(x) matrix element receives 
contributions from T-even PDFs, while the <&g(x) matrix element receives contri- 
butions from T-odd PDFs, see Ref. [14] for a detailed study of this. For the single 
weighted results, thus, the behavior under time-reversal can be used to identify the 
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process dependent parts and we find 

f -L(i)[c/] r , 



(x) — c\j ' f 1T ■* (x) , 



(18) 



J IT 



where transverse weighting for PDFs is defined as 



(19) 



for weighting with n transverse moments. The fact that de pen ding on the process 
are the correlators to be used, leads to the sign change for the T-odd func- 
tions in such processes or to more complex factors if more complex gauge links are 
involved ^1 The importance of Eq. (|18[) is the appearance of a universal function 
with calculable process (link) dependent numbers rather than many process depen- 
dent functions that are somehow related. For gluon TMDs, there are already for 
single weighting two functions Tq' 1 ' and hence two different gluonic pole factors 
Cjj ^' , because there are two ways to construct color singlets from the (in that 
case) three gluon fields that arc involved using the f a b c or d a b c structure constants. 
The appearance of two different gluon Sivers functions was poi nte d out in Ref. 1171 

1 1 Ql m 

The situation with universality for fragmentation functions is easier because 

11 QI90I91 1 

the gluonic pole matrix elements vanish in that case 1 1 1 I Nevertheless, there 
exist T-odd fragmentation functions, but their QCD operator structure is T-even. 
These T-odd functions then appear in the parametrization of Hence, there is 
no process dependence from gluonic pole factors. 

The use of transverse moments in the description of azimuthal asymmetries via 
transverse momentum weighting of the cross section can be extended to higher 
moments involving higher harmonics such as cos(2y>). Also here process dependence 
may come in from double gluonic pole matrix elements Q'qq, which are twist four 
operators. This affects studies that involve the quark TMD h^(x,p T ) (pretzelocity 
distribution) in Eq. (fTTj) or the gluon Boer- Mulders function h 1 9 (x,p T ) (linearly 
polarized gluons in unpolarized targets) in Eq. (|13j). For instance, for quarks one 
finds for the simplest gauge links, 



more complex gauge links one actually needs a third (universal) pretzelocity second 
transverse moment because there is another possible color structure 

4. TMDs of definite rank 

An interesting possibility to obtain universal TMDs is to start with a parametriza- 
tion that involves the symmetric traceless tensors p" 1 "'"™ of rank m, such as 




Pti Pt PtPt ^Pt 9t 



(21) 
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Depending on the rank different correlators come in, involving operator combina- 
tions of gluons, covariant derivatives and ^4-fields. Minimizing the twist we have 

^(x,p T ) = *M)+»C^^^ > ^)+^C^ ie ^4g 0ie ( a ; > 4) + ... 

+ ^ + ••■. (22) 

with a summation over the color structures c. The reproduction of the transverse 
moments provides the proper identification of universal TMD functions, 



${x,pl) = <yfi(x,p 2 T ) + S L gi{x,pl)~f 5 + h 1 (x,pl)^ 5 $ T j — , (23) 
^$M^^) = {^L(-,^)^ 2 ^-3iT(^^)^7 5 }f (24) 
^ $U*,pI) = l ~{- £P/ ^f rCT + ihi(x, P l) |}|, (25) 



3 



^*Ux,pl) = h£ A \x,£) PTij ^ (26) 

P ^% G ,{x,pl) = ±htr\^l) |, (27) 

PlU^j (xv 2 )- — h ±{B2) (x v 2) PliSj^Z t (28) 

^*^(*,l£) = 0. (29) 

We note that the rank zero functions in Eq. (|23[) depend on x and p2, and in- 
volve traces, to be precise gi(x,p 2 ) = gf^(x,p 2 ) and hi(x,p 2 ) — h^(x,p 2 ) — 
(p 2 T /2M 2 )h^ U] (x,p 2 T ). As remarked before, for the pretzelocity there are three uni- 
versal functions with in general 

h^ u \x,pl) = h^ A \x,p 2 T ) + Cfil h$ m \x,p 2 T ) + Cf^ 2 h^ B2 \x,p 2 T ). (30) 

For the simplest gauge links we have Cqq 1 = 1 and Cqq 2 = 0, which shows 

e.g. that hf^ IDIS \x,p 2 ) — hf^ Y \x,p 2 ), but that for other processes (with more 
complicated gauge links) other combinations of the three possible pretzelocity func- 
tions occur. For a spin 1/2 target the abov e set o f TMDs is complete. There are 
no higher rank functi ons. For a spin 1 target ESiLMI an( j f or gi uonSj there are higher 
rank functions | 22 | 25 | p or fragmentation correlator there is for rank 2 only a 
single (T-even) pretzelocity function H^ T {z : k 2 ) appearing in the parametrization 
of the correlator Aqq(x,Pt). The rank of the various correlators is shown in Table[T] 
with the results for nucleon TMD PDFs summarized in Tables [2] (unpolarized) and 
El (polarized) and those for nucleon TMD PFFs in Tables U and El 
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Table 1. The correlators of definite rank in the full TMD correlator, ordered in columns 
according to the number of gluonic poles (G) and ordered in rows according to the number of 
contributing partial derivatives (d = D — A). The rank of these operators is equal to the sum 
of these numbers. Their twist is equal to the rank + 2. 



GLUONIC POLE RANK 





1 


2 


3 








^ C GGG,c ®GGG,c 


<5>g 




* 2 C l GG,c * {dGG h c 





















Table 2. TMD PDFs for an unpolarized or 
spin target assigned to the quark correlators 
as given in Table [T] 



PDFs FOR SPIN HADRONS 


h 



















Table 3. Assignments of TMD PDFs for a po- 
larized spin 1/2 target. 



PDFs FOR SPIN 1/2 HADRONS 


91, hi 


/it 




9it, h^ L 







"it 



Table 4. Assignments of TMD PFFs for un- 
polarized or spin hadrons. 



PFFs FOR SPIN HADRONS 


Di 






Hi 













Table 5. Assignments of TMD PFFs for po- 
larized spin 1/2 hadrons. 



PFFs FOR SPIN 1/2 HADRONS 


Gi, Hi 






Git, H 1L , D^ T 






H 1T 







5. Conclusions 

We have introduced quark TMD correlators of definite rank in Eq. (|22| . In this new 
decomposition, we have made an expansion of the quark correlator into irreducible 
tensors multiplying correlators containing operator combinations of gluons, covari- 
ant derivatives and A-ficlds, the latter in the combination id = iD — gA. In the 
decomposition gluonic pole factors contain the gauge link dependence, which are 
calculated from the transverse moments. These factors also give the process depen- 
dence, which is determined by the gauge link structure. The correlators of definite 
rank, in turn are parameterized in terms of the universal TMD PDFs depending on 
x and pj, such as given by Eqs. (|23p -(|29 p . The process dependence for a particu- 
lar TMD PDF is in the same gluonic pole factors that appear in the expansion in 
Eq. 422]). 

An analysis for a quark spin 1/2 target shows that the process dependence 
is not strictly confined to the T-odd functions, such as the Sivers or the Boer- 
Mulders functions. In fact, we have shown the existence of three T-even pretzelocity 
functions. For fragmentation the TMD PFFs are already universal since gluonic pole 
matrix elements vanish for fragmentation correlators. Future work will be focused 
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on the study of universality for higher spin targets and gluon TMD PDFs l ^ l ^ l 
While for a spin 1/2 target one has at most rank two TMDs, one has for higher 
spins and gluon TMDs also higher rank functions, while also the color and gauge 
link structure is richer. 
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